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Abstract 

A geometrical characterization of robustly separable (that is, re- 
maining separable under sumciently small variiations) mixed states 
of a bipartite quantum system is given. It is shown that the density 
matrix of any such state can be represented as a normal vector to a 
hypersurface in the Euclidean space of all self-adjoint operators in the 
state space of the whole system. The expression for this hypersurface 
is provided. 

1 Introduction 

Entanglement turned out to be a crucial resource for quantum computation. 
It plays a central ròle in quantum communication and quantum computation. 
A considerable effort is being put into quantifying quantum entanglement. 

It seems natural to focus the effort s on quantifying entanglement itself, 
that is, describing the impossibiliti/ to prepare a state by means of LOCC (lo- 
cal operations and classical Communications). One may, although, go another 
way around and try to quantify separability rather than entanglement: this 
turned out to be applicable for building combinatorial entanglement patterns 
for multipartite quantum systems [3]. 

In this paper I dwell on the case of bipartite quantum systems. A state 
of such system is called separable if it can be prepared by LOCC. In terms 
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of density matrices that means that p, its density matrix, can be represented 
as a mixture of pure product states. According to Carathéodory theorem, 
the number of this states can be reduced to n A where n is the dimension of 
the state of a single particle. 

The idea to replace finite sums of projectors by continuous distributions 
on the set of unit vectors is put forward making it possible to provide a 
geometrical characterization of separable mixed states of a bipartite quantum 
system. To consistently describe the result presented in this paper recali some 
necessary definitions. 

Bàsics. A density matrix p in the product space 23 is called factorizable 
if it is a tensor product of density matrices, p = p <E> p' ■ If p is a convex 
combination of factorizable operators, it is said to be separable 

a 

A crucial feature of quantum mechanics, the phenomenon of quantum en- 
tanglement, stems from the fact that there exist density operators in the 
product space which are NOT separable, they are called entangled. A 
density operator p is called robustly separable if it has a neighborhood 
U in C such that all operators p' G U are separable. 

A brief account. In the Euclidean space £ of self-adjoint operators acting 
in the tensor product space 23 = Ti <g> Ti' we define a real-valued, positive 
functional K, : £ — > M + as follows 

K(X) = JJe^ x ^dS n dS' n 

where the integration is taken over the torus — the Cartesian product of unit 
spheres in Ti, Ti', respectively, and consider the hypersurface /C C £ 

/C = {X g £ | K(X) = 1} 

Then 

• all robustly separable density operators in Ti are in 1-1 correspondence 
with the points of /C 
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• the density matrix associated with a point X e K, is the normal vector 
to K, at point X. 

2 Continuous optimal ensembles 

To make the account self-consistent, begin with necessary definitions. A den- 
sity operator is a non-negative self-adjoint operator whose trace equals to 
1. In particular, for any unit vector \<fi) the one-dimensional projector \<j))(<t>\ 
is a density matrix. Note that for any set of density operators p a the convex 
combination J2 a P» * s a l wa Y s a density operator. 

The set of all self-adjoint operators in TC = C n has a natural structure of 
a real space M 2n , in which the set of all density matrices is a hypersurface, 
which is the zero surface T = of the affine functional T = Tr X — 1 . 

In this paper a geometrical characterization of separable bipartite den- 
sity operators is provided. It is based on the notion of continuous ensembles. 
Generalizing the fact that any convex combination of density operators is 
again a density operator, we represent density operators as probability dis- 
tributions on the unit sphere in the state space TC of the system. Let us pass 
to a more detailed account of this issue beginning with the case of a single 
quantum system. 

Let TC = C n be a n-dimensional Hermitian space, let p be a density matrix 
in TC. We would like to represent the state whose density operator is p by an 
ensemble of pure states. We would like this ensemble to be continuous with 
the probability density expressed by a function p((fi) where (fi ranges over all 
unit vectors in TC. 

Technical remark. Pure states form a projective space rather than the 
unit sphere in TC. On the other hand, one may integrate over any probabilistic 
space. Usually distributions of pure states over the spectrum of observables 
are studied, sometimes probability distributions on the projective spaces are 
considered 0. In this paper for technical reasons I prefer to represent en- 
sembles of pure states by measures on unit vectors in TC. I use the Umegaki 
measure on CB n — the uniform measure with respect to the action of U (n) 
normalized so that J CB dS n = 1. 
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2.1 Effective definit ion 



The density operator of a continuous ensemble associated with the measure 
on the set CB n of unit vectors in 7ï is calculated as the following 
(matrix) integral 



p = J n{cj>) |0X0| dS n (2) 

where is the projector onto the vector (0| and dS n is the above men- 

tioned normalized measure on CEL: 



dS n = 1 (3) 

<t>&<CB n 

Effectively, the operator integral p in (J2J) can be calculated by its matrix 
elements. In any fixed basis {|ej)} in 7í, each its matrix element pij = 
(ej| p \ej) is the following numerical integral: 

pij = (ei\p\ej) = J p,{4>) (e i |0)(0|e i ) dS n (4) 

<t>eCB n 

2.2 Optimal ensembles 

We need to solve the following variational problem. Given a functional Q on 
L 1 (CS n ) and given a density matrix p in H, find the distribution p on the 
set CB n of unit vectors in Ti such that 

/ M</0 10X01 rfS « = p 

(5) 

<5(/i) — * extr 
We shall consider functionals Q of the form 

Q(p) = j q(p(<P))dS n (6) 

(t>&CB n 

then, according to (j4*|) . the variational problem (JHJ) reads 
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/ n((f>) (ei | 0) (0 | ej) dS n = p iS 
/ q{n{4>)) dS n -> extr 

<j)GCB n 

Solving this variational problem by introducing Lagrangian múltiples Xy we 

get 

- Çx^e, 10x01è,-) = o (?) 

Combining the Lagrange múltiples into the operator X = Ylij^ij l e j)( e íl 
turns the equation (J7J) to ç'(/i(0)) = (0|X|0). Then, denoting by / the 
inverse of q' we write ((Zj) as 

M0) = /((0|X|0)) (8) 
and the problem reduces to finding p from the condition 



M0) 10X01 dS n = p (9) 
4>eCB n 

which according to (JHI) and (jlj) can be written as 

fel p|e,-) = f f ((0| X |0)) |0)(0| dS n (10) 



It follows from (J2J) that the coefficients can be chosen so that Xq. = X^. 
That means that the problem of finding the optimal ensemble reduces to 
that of finding the coefficients of a self-adjoint operator, that is, to finding 
n 2 numbers from n 2 equations. 



2.3 Geometrical interpretation 

The equation (|10|) can be given a direct geometrical meaning. Let C ~ R n2 be 
the space of all self-adjoint operators in 7í. Let / : R — ► R be a differentiable 
function. Consider the real valued functional F : £ —>■ R defined as 

F(X) = J /((0|X|0)) dS n (11) 

<f>€CB n 
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which is well-defined as the set CB n is compact. Fix a basis {e^} in TC, 
then any X e C is defined by its matrix elements = (ej|X|efc), so 
(0| X |0) = J2ik (0 | 6j) (efe | 0). Then the expression (fTT]) can be treated 
as an integral depending on the set of parameters {Xa-}. We may consider 
the derivatives of F(X) with respect to these variables, calculate them 

d f d 



8xj^= I 4- t ( fmxm ) ds " = 

<t>eCB n 

f ((0|X|0))(0|e,,) (e fe |0) dS n = (12) 



4>eCB n 

= (e fc | | f ((0| X |0)) |0)(0| dS n | ei ) 
<t>eCB n 

So, the gradient of the functional F is the operator which can be symbolically 
written as 



VF = J /' «0| X |0» Í0X0| dS n (13) 
<t>eCB n 

and effectively calculated using (JT2J). 

2.4 Optimal entropy ensembles 

Let us specify the form of the optimality functional in (JHJ) assuming it to be 
the differential entropy of the appropriate distribution: 

g(/i) = —/i ln \i (14) 
then q' — — (1 + hi fi) and we have the following / for (jlOj) 

f(x) = e ~ {1+x) 

Introduce, as in (|llj). the functional /C : C — > R on the set of all self-adjoint 
operators in 7i (the minus sign and the unit summand are omitted here being 
a matter of renormalization): 

K(X) = J e^^dSn (15) 



6 



Note that p(X) = J e ^' x '^ \4>)(<j>\ dS n is always a positive operator, 

<j>£CB n 

then 



K{X) = Tr J e W x W \<f))((j>\ dS r 



<j}€CB n 

is a condition which defines a full-range density matrix p(X) in 7ï. On the 
other hand, the condition K{X) = 1 defines a hypersurface in the Euclidean 
space C. Together with the fact that (e x ) = e x and ()13|) we come to the 
following 

Statement. Any full-range density matrix in p is associated with a point 
on the hypersurface fC(X) = 1 and the entries of p are calculated as the 
components of the gradient: 



p = V/C (16) 



2.5 The existence 



Why optimal entropy ensembles do exist for all full-range density matrices? 
First note that for any full-range density matrix p = YlPk l e fc)( e fc| there 
are infinitely many continuous ensembles (=probability measures on CB n 
in our setting) associated with it. An example of such distribution is p = 
Pk |e fc ){e fc | = / fi(<p) dS n with 

as it follows from jïïj. Here L is a parameter, such that L > po ^+i) where 
p > is the smallest eigenvalue of p. Any probabilistic density p whose 
support is CB n is a point in the interior of the simplex of all probabilistic 
measures on <CB n . For each probabilistic measure on CB n its differential 
entropy can be calculated. The differential entropy is, in turn, a concave 
function in the affine space of probability distributions. Therefore if we have 
an affine subset of of probability measure on CB n , the differential entropy 
takes its maximal value in the interior of the simplex of probability measures. 
Now return to the condition in (JHJ) — we see that it is affine. Therefore, if we 
know that there exist at least one continuous ensemble representing p (but 
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we know that as mentioned above), that means that there exist a maximal 
entropy ensemble representing X, hence it has the representation (JTBJ). 

3 Bipartite systems 

Consider two finite-dimensional quantum systems whose state spaces are 
TC, TC'. The state space of the composite system is the tensor product 03 = 
TC <S> TC'. Denote by £ = C ® C the space of all self-adjoint operators in 03. 

3.1 Continuous ensembles in bipartite case 

Let p be a robustly separable density matrix in the product space TC ® 
TC'. Then it can be represented (in infmitely many ways) as a continuous 
ensemble of pure product states. Carrying out exactly the same reasoning 
as in section 1231 we conclude that among those continuous ensembles there 
exists one having the least differential entropy, this will be the ensemble 
we are interested in. Like in section EU formulate the variational problem. 
Let p be a density operator in a tensor product space 03 = TC <S> TC'. The 
task is to find a probability density /i(0</>') defined on the Cartesian product 
T = CB n x CB n of the unit spheres in TC,TC', respectively. 



I Q(/i) — > extr 

Proceeding exactly in the same way as with single particle, we get the 
following representation: 



for some self-adjoint operator X in C whose existence is guaranteed by the 
same reasons as in section 12.51 Why such X does not exist for entangled 
density operators? The reason is that the set of probability distributions 
among which e^^ x ^' is optimal is simply void in the entangled case. 




(18) 




(19) 
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3.2 Geometrical characterization of robustly separable 
quant um states 

Now we pass to the main result of this paper. Suppose we deal with a tensor 
product of two Hilbert spaces 7í, Tt', each of dimension n. Consider the space 
C of all self-adjoint linear operators in the tensor product 23 = H®H', being 
a Euclidean space of dimension n 4 . For any X G C we can always calculate 
the integral 

K(X) = J e^'W^dSndS'n (20) 

which is always well-defined (as an integral of a bounded function over a 
compact set), positive (as the exponent is always positive) functional from 
C to R+. Consider the hypersurface K in C defined by the equation 

K = {X eC\IC(X) = 1} 

In any point of C the gradient V/C can be calculated. In particular, at 
any point X of /C the gradient V/C will be a normal vector to /C. The 
surface /C is something given once and forever, it depends only on the di- 
mensionality of the state space. For any X such that JC(X) = 1, we can 
calculate the gradient p(X) = V/C \x at point X Fix bases {ej}, {e-,}, then 
X = J2u> kk' Xü'kk' \ii')(kk'\ and the expression (fTïïj) for the operator p has 
the following form: 

Pwkk' = V/C = — — (21) 

OAn' kk' 

Conversely, given a robustly separable bipartite density matrix p, we 
know that it can be represented as a convex combination of product states: 
p = ^2p a p a ® p' a - Each p a can be, in turn, represented as a non-vanishing 
probability distribution (|17|) . Then exactly the same reasoning as in section 
12.51 can be carried out and there is a point X on the surface /C associated 
with p. So, together with (|2*Tj) . we have the main result: 

{robustly separable states} <-> {the points of /C} (22) 
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Summary 



A geometrical interpretation of robustly separable density operators of a 
bipartite quantum system with the state space 23 = 7í®7í' is provided. They 
are represented as normal vectors to the hypersurface K, in the (Euclidean) 
space £ of self-adjoint operators in 23 defined by the following equation: 



where the integration is performed over the set of all unit product vectors 
(4>4>'\ € 23. Each point X 6 K, is a self-adjoint operator, the parameter of 
the probability distribution on the set of unit vectors which gives a density 
operator p. Furthermore, the normal vector to /C at point X is p itself: 



The final remark. Given a density matrix p in 23, a question arises if it is 
separable or not. When the dimension of at least one of spaces TC, W is 2, this 
question was given an effective answer — the positive partial transpose (PPT) 
criterion due to Peres- Horodecki was suggested [T]. The criterion states that 
p is separable if and only if its partial transpose p T2 remains non-negative 
matrix. In higher dimensions PPT is only a necessary condition for a state 
to be factorizable as there exist entangled density matrices whose partial 
transpose if positive. 

Although a geometrical characterization of robustly separable density ma- 
trices is provided, it does not solve (directly, at least) the 'inverse problem'. 
Nevertheless, the continuous ensemble method presented in this paper seems 
to be helpful for tackling the inverse problem as well. This issue is addressed 
in the next paper on continuous ensembles. 
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